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Abstract 

We employ the Landau-Lifshitz-Bloch (LLB) equation to investigate current-induced domain 
wall motion at finite temperatures by numerical micromagnetic simulations. We extend the LLB 
equation with spin torque terms that account for the effect of spin-polarized currents and we find 
that the velocities depend strongly on the interplay between adiabatic and non-adiabatic spin 
torque terms. As a function of temperature, we find non-monotonous behavior, which might be 
useful to determine the relative strengths of the spin torque terms experimentally. 

PACS numbers: 72.25.Ba Spin polarized transport in metals. 75.10.Hk Classical spin models, 75.60.Ch 
Domain walls and domain structure 
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I. INTRODUCTION 



Magnetic nanostructures in external magnetic fields as well as under the influence of 
spin-polarized currents have become interesting research fields in recent years due to funda- 
mental novel effects that occur for geometrically confined spin structures, such as domain 
walls.- Current-induced domain wall motion has been suggested as an alternative to the 
use of external magnetic fields to induce switching, opening the possibility of simple device 
fabrication making field-generating strip lines redundant. While current-induced domain 
wall motion is experimentally well established,- 1 ^ the underlying physical mechanisms are 
not completely understood yet and in particular the importance of the adiabatic and the 
non-adiabatic spin torque terms as well as domain wall transformations for high current den- 
sities are highly debated.-^ Furthermore, the influence of temperature on the effects has 
so far been neglected in the K calculations and so the experimentally found temperature 
dependence of for instance the critical current densities is so far not understood.- 

To theoretically predict the behavior of a spin texture under current, one can numerically 
solve the Landau-Lifshitz-Gilbert (LLG) equation and computer simulations can be per- 
formed using either a micromagnetic model or a classical atomistic spin model. Spin torque 
effects can be taken into account by including the adiabatic and the non-adiabatic torque 
terms.-^ 1 ^ Due to the computational expense of atomistic simulations, system sizes are 
restricted to a nanometer range, so that micromagnetic approaches are desirable. However, 
conventional micromagnetic calculations for larger system sizes lack the correct description 
of temperature effects because of the assumption of a constant magnetization length. An 
alternative approach that has only recently started to be used widely to investigate realistic 
systems sizes including temperature effects is to employ the so-called Landau-Lifshitz-Bloch 
equation.— This equation forms the basis for micromagnetic calculations at elevated tem- 
peratures using a macro-spin model where longitudinal relaxation processes are taken into 
account^ but so far the LLB equation has only been studied without the spin torque terms. 

In this paper, we extend the LLB equation of motion by adding the spin torque terms and 
we study domain wall motion under the influence of current and at variable temperatures. 
We determine the domain wall velocities and find that they exhibit a strong dependence on 
the temperature. Furthermore, by the interplay between the adiabatic and the non-adiabatic 
spin torque the resulting onset of domain wall transformations (Walker breakdown) is very 
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sensitive to the temperature. 



II. MODEL 

A. Landau-Lifshitz-Bloch equation 

While in the LLG equation at K the length of the macro-spins stays constant, for finite 
temperatures an equation of motion for macro-spins allowing for longitudinal relaxation 
was derived by Garanin^ within mean-field approximation from the classical Fokker-Planck 
equation for atomistic spins interacting with a heat bath. The resulting "Landau-Lifshitz- 
Bloch equation" has been shown to be able to describe linear domain walls, a domain wall 
type with non-constant magnetization length.— Furthermore, the predictions for the 
longitudinal and transverse relaxation times have been successfully compared with atomistic 
simulations^ as well as rapid heating experiments.— Therefore, we now employ this equation 
to study the thermodynamics as well as the excitations of macro-spins due to currents. 

The LLB equation can be written in the form 



where rhi is a to its zero temperature value normalized spin polarization and 7 the gy- 
romagnetic ratio. The magnetization is not assumed to be of constant length and even 
its equilibrium value, m e , is temperature dependent. Hence, besides the usual precession 
and relaxation terms, the LLB equation contains another term which controls longitudinal 
relaxation. 

The LLB equation is valid for finite temperatures and even above the Curie temperature 
Tc though the damping parameters and effective fields are different below and above To 
an and a± are dimensionless longitudinal and transverse damping parameters. For T <Tq 
they are an = 2AT/(3Tc) and aj_ = A(l — T/(3Tq)). For T >Tq the damping parameters 
are equal, a± = a\\ = 2AT/(3Tc). Here, A is a microscopic damping parameter which 
characterizes the coupling of the individual, atomistic spins to the heat bath. In the limit 
T — > the longitudinal damping parameter an vanishes and with a± = A the LLB equation 
evolves into the usual Landau-Lifshitz (LL) equation. 




(1) 
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The effective fields of the LLB equation are the derivative = — j^oji- of the free 
energy density /. The total field is given by^ 
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■ 1- ^ \rfii T<T C 
Hl s = Hl + Hl x +e x \\ ' ~ (2) 

H( i+ i(S)h T ^ T c 



with the biaxial anisotropy field, 

HI = (hnfa + m %) > (3) 

which makes the x-axis the easy axis, the y-axis the intermediate axis and the z-axis the 
hard axis of the model. The exchange field is 

were A is the lateral size of the discretized cells, A is the temperature dependent exchange 
stiffness, and M s ° is the zero temperature saturation magnetization. The susceptibilities xi 
are defined by xi = dmi/dBi with / =||, _L. Note, that at low temperatures the perpendicular 
susceptibility x± is related to the temperature dependent anisotropy constant K via x± = 
M®m%/(2K).— We use functions for Xi(T), m c (T), and A(T) as calculated before for the 
spin model (for details seeii^) but rescaled to reflect a ferromagnetic material with a Curie 
temperature of 1043 K and an M s ° of 10 6 A/m. Furthermore, we normalize the perpendicular 
susceptibility such that its value at OK, x_i_(T = 0) = M s °/2i^ x corresponds to an anisotropy 
of K x = 10 5 J/m 3 and 2K y = K x . These functions are shown in Figs. H]and[2j 

B. Spin torque in the Landau-Lifshitz-Bloch equation 

Throughout this paper we will consider a one-dimensional model of a domain wall. An 
established approach for the effect of a spin-polarized current in the x-direction on a domain 
wall is presented in Ref. (]. In these studies the interaction between electron spins and 
magnetization has been treated by additional spin torque terms 

dS -> dS 
T = -u x — + pSxu x —, 5 
ox ox 

where S is a unit vector representing the direction of the magnetization. The first contri- 
bution to the spin torque is called the adiabatic term. It can be derived from an additional 
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FIG. 1: (Color online) Equilibrium parallel and transverse susceptibility vs. temperature deter- 
mined as explained in the text. The inset shows the transverse susceptibility vs. temperature on 
a larger scale. 
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FIG. 2: (Color online) Exchange stiffness A and reduced equilibrium magnetization m e vs. tem- 
perature determined as explained in the text. 

term in the magnetic free energy that takes into account the coupling of the magnetization 
to the spins of the electrical current,— representing an adiabatic transfer of angular mo- 
mentum to the magnetization. In the adiabatic limit the spin polarization of the current is 
always oriented along the local direction of the magnetization. The second contribution is 
the non-adiabatic term that reflects the mistracking of the direction of the conduction elec- 
tron spins with respect to the magnetization. It appears to play a role similar to the Gilbert 
damping term. Recent micromagnetic numerical investigations^'&^^i using a modified 
LLG equation including these terms have given a qualitative insight into the roles played by 
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these two torque terms. 

In the case of the current flowing in the re-direction, the magnitude of the effective spin 
current u x is given by w x = j x /M where M is the magnetization and j x the spin current 
density j x = ^Pje/e which is proportional to the electrical current density j e and to 
the polarization P. Here, /xb is the Bohr magneton and e the magnitude of the electron 
charge. Normally the spin current density is assumed to be temperature independent. In 
the following, we extend the model so that a temperature dependent spin current j x (T) is 
taken into account. Under the assumption that the spin polarization is proportional to the 
magnetization, P(M) = P°m, and with M = M s °m, u x is u x = - '^fj^ for all temperatures. 

Under these assumptions the spin torque terms can be expressed in terms of the reduced 
magnetization m as 

- dm 3 _ dm 

T = -u x — H m x « x — . (6) 

ox m ox 

Note, that the spin torque now is temperature dependent via the variable m which, within 
the framework of the LLB equation, is no longer a unit vector of constant length. For the 
same reason, the adiabatic torque term can no longer be expressed as a double cross product. 
Instead it is 

dm u y / _ , ^ dm s , _ dm , 



which means that the adiabatic term, within the LLB equation, gives rise to an additional 
longitudinal spin torque term which vanishes in the LLG equation due to the assumption of a 



constant length of the magnetization vector. As pointed out in Ref. |21| this term corresponds 
to the spin accumulation and in the metal systems considered here, it constitutes usually a 
small effect. 

It is not yet fully clear whether the spin torque term T should be added to the Landau- 



Lifshitz or the Landau-Lifshitz-Gi 
see the discussion in Refs. 
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bert form of the equation of motion (for more details 
23IJ24J ) . The same problem arises with the LLB equation. In 
the following, we extend the LLB equation with both forms of the damping, the one after 
Landau and Lifshitz as well as the one after Gilbert. The LL form of the LLB equation 
(Eq. [T]) is the original one as derived by Garanin.— Eq. [T] now reads with the additional spin 
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torque terms from Eq. [6] 

m i= - ry m . x ff eS _ m . x irriiX H c{{ 1 

m i \ J 

m\ \ J 

dm /3 hL _ drhi 

- u x — 1 mi x u x ^— . (8) 

ox mi ox 

Neglecting terms of the order of a 2 the LLB equation can be transformed into an equation 
with a damping term following Gilbert. Adding the same spin torque terms T to this form 
of the LLB equation yields 

rhi — — •yrhi x H l cS H T rhi x m, 

m i 



+ 2" m * ' -"cflF ) m 'i 

m\ V / 



- M x — 1 TO; X M x — — . (9) 

ax rrti ax 

Note, that we use the notation /?g, for the non-adiabatic prefactor only for convenience 
in order to distinguish in the following between the LL and LLG form of the LLB equation. 

In the next step we transform Eq. [9] into an explicit form so that we are able to compare it 
with Eq. [8] and also since an explicit equation is more convenient for a numerical treatment. 
This explicit equation can be derived once again neglecting terms of the order a 2 and a/3 
and it is given by 



rhi = - irfii x H! s r mj x [rhi x H, 

m\ V 

+ ^(rhi-K^m 



drhi (fa a±\ _ drhi 

u x — h ^ m,x« x — . (10) 

ox Vm,- mfs ox 



The only difference between Eq. [8] assuming Landau-Lifshitz damping, and Eq. [10] assuming 
Gilbert damping, is the prefactor of the last term. Eqs. [8] and [10] are mathematically identical 
for — — oc±/rrii. Note, that at OK both equations evolve into the well established 
explicit versions of the LL respectively LLG equation with spin torque terms.— i22a2i 

III. ANALYTICAL MODEL 

In the following, we present analytical calculations for a one-dimensional system aligned 
along the x-direction. We extended well-established, analytical OK calculations^ 1 ^ to ele- 
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vated temperatures using the LLB equation as equation of motion. Note, that within the 
LLB approach domain wall profiles change with increasing temperature first from circular 
to elliptical and later on to linear.— ii&M However, these effects occur only at higher temper- 
atures close to Tq (depending on the strength of the anisotropy), and are neglected in our 
analytical calculations. Hence, we assume a temperature independent domain wall type, 
which here is a transverse domain wall along the x-direction, which does not change its 
spin structure dynamically. This assumption, is later tested by comparison with numerical 
calculations without these approximations. 

The assumed domain wall profile is as usuat^ 1 ^ 



m = —m. 



oosh(f) v cosh(f) 
with the temperature dependent domain wall width 




^ T ^ = —o\l ml{T) (1 + sin20) ( 12 ) 
and the out-of-plane angle 0. In the low temperature limit this equation has the well known 



form, A = ^2A/{K{1 + sin 2 0)). 

The equations of motion are calculated as described in Ref-O. We assume Gilbert damping 
first, simply because it is the more common assumption in connection with spin torque 
calculations. For the domain wall profile above the equations of motion for the position of 
the domain wall Xq and the angle are: 

m c A / -f 

%n = — - — sin20 + M x (13) 
sin 20 — (14) 



4x± m e A 

These equations are calculated based on the original procedure (see^ and references therein). 
The corresponding Walker threshold u}^^— can be calculated under the assumptions that 
= and that the domain wall width A(0, T) reaches its minimum A min at = ±7r/4. 
The Walker threshold is then given by 

u G - —^—A a±me (15) 
WaIker ~4x± Ha±-/3 G m e | 

and the average domain wall velocity is 

P G u x m e m e A min 7 



(«>( 



m e A min 7 f, i£ ~ 
± 4 y, ) _1 ' (16) 

A-L V "Walker 




FIG. 3: (Color online) Walker threshold ^walker according to Eq. [15] vs. the reduced temperature 
T/Tc for different values of the non-adiabatic prefactor (3q as well as 

with the plus sign for (a± — (3Qm e )u x > and the minus sign otherwise. This equation 
contains a contribution which is linear in the current (in the non-adiabatic case) and a 
second square-root contribution above the Walker threshold. The temperature dependence 
is included in the temperature dependent parameters A(T) , m e (T) , x±{T) , and a±(T). At 
zero temperature these results are identical to those gained with the LLG equation.-*^ The 
longitudinal spin torque does not affect the analytical results since we assume a constant 
domain wall type. Very close to the Curie temperature this assumption is no longer valid 
and deviations can be expected, which are beyond the scope of the current investigation. 

In Fig. [3] the Walker threshold according to Eq. [15] is shown as a function of the reduced 
temperature (T/Tc). In all cases the Walker threshold vanishes at the Curie temperature. 
In the adiabatic case (/3q = 0) and in general for (3q < A = a±(T = 0) it decreases with 
increasing temperatures. For (3q = A = 0.02 the Walker threshold diverges approaching zero 
temperature. The behavior for (3q > A is similar, but here the Walker threshold diverges at 
a finite temperature. As we will see in the following, the understanding of the temperature 
dependence of the Walker threshold is the key for understanding the temperature dependence 
of the domain wall velocity. 

The equation of motion for Landau-Lifshitz damping (Eq. [TOT) can be calculated and solved 
in the same way as for the Gilbert damping with corresponding results for the Walker thresh- 
old and the average domain wall velocity for = ~ a±/m c . Note, however, that since 
this transformation is temperature dependent, the temperature dependence of both, Walker 
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threshold and domain wall velocity are different for Landau-Lifshitz and Gilbert damping, re- 
spectively. The temperature dependence of the Walker threshold assuming Landau-Lifshitz 
damping is also shown in Fig. [3J The main difference is that approaching the Curie temper- 
ature the Walker threshold does not vanish but diverges. In the limit of low temperatures, 
however, the Walker threshold converges to the one following Gilbert damping. 

IV. NUMERICAL MODEL 

By means of computer simulations, temperature dependent domain wall velocities were 
calculated for a one-dimensional system of 512 nm length, discretized with one nm cell size. 
The initial magnetization configuration was a planar domain wall positioned in the middle 
of the chain with the temperature dependent profile and width given by Eqs. [11] and [12] 
with = 0. For OK the domain wall width varies between A max = 20 nm (0 = 0) and 
Amin = 14.26 nm (0 = 7r/4). At the ends of the system the spins were fixed as boundary 
conditions in the x-direction and -x-direction, respectively. To minimize the influence of 
these boundary condition the domain wall was only allowed to move within a range of 60 nm 
from the center of the system. When the domain wall moved out of this interval, it was 
shifted back along the x-coordinate and re-positioned at the opposite side of the interval. 
The domain wall velocity was calculated from the derivative of the spatially averaged x- 
component of the magnetization versus time. The numerical time integration of Eqs. [H]and 
[TUJwas carried out using a Heun-method^^ with 1.8 fs time step size. 

V. RESULTS 

A. Adiabatic spin torque effect 

In the following, first the pure adiabatic spin torque effect (/3q = 0) is discussed. 
For zero temperature, this effect was already investigated previously based on the LLG 
equation.— ii£i2&2£ It was found that for low effective spin currents u x the domain wall moves 
along the wire until it reaches a maximum displacement where it stops. At the same time, 
the magnetization of the domain wall is tilted out of the easy plane up to a maximum out- 
of-plane angle. This behavior can be explained by an analysis of the different terms of the 
extended LLG equation (Eq. [TOl) : the first spin torque term which moves the domain wall 
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FIG. 4: (Color online) The top figure shows the average domain wall velocity (v)g calculated from 
the LLB equation with Gilbert damping vs. the effective spin current u x for different temperatures. 
The solid line for T = K represents the case where the LLB and the LLG equation are identical. 
The bottom figure shows {v)q vs. the reduced temperature T/Tq. The lines represent the analytic 
solution (Eq. [To]) and the points are from numerical simulations. The calculations are for the 
adiabatic case, (3q = and A = 0.02. 



is balanced by an "internal" torque due to the anisotropy contribution to the effective field. 
The displacement of the domain wall in x-direction is eventually stopped by the precessional 
term acting in the opposite direction, while the second spin torque term which is responsible 
for tilting the magnetization out of the easy plane is balanced by the damping term. 

It was even analytically predictecr^^ that the averaged domain wall velocity as a 
function of the effective spin current remains zero unless the current exceeds a critical value 
u c . This predicted critical current was also found in atomistic simulations at K.— Below the 
critical effective spin current (u x < u c ) no continuous domain wall motion is observed while 
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above the critical current the spin torque term can no longer be balanced by the anisotropy. 
Consequently domain wall motion occurs in addition to a precession of the magnetization 
around the x-axis. 

We find the same behavior in the extended LLB equation. This can be seen in Fig. H] 
where the averaged domain wall velocity is shown as a function of the effective spin current 
u x for different reduced temperatures T/Tq. However, it is found that the critical effective 
spin current is temperature dependent, following the equation 

U c = — —&mm (17) 

(Eq. [15] for (3q = 0). For larger temperatures the critical current decreases since thermo- 
dynamically the anisotropy decreases. This leads to the fact that domain walls at higher 
temperatures are faster than at low temperatures. 

Furthermore, in the bottom part of Fig. 0] the averaged domain wall velocity is shown 
as a function of the reduced temperature for different values of the effective spin current. 
In the limit of low current a critical temperature T* can be identified. For T < T* no 
continuous domain wall motion is observed while for T > T* domain wall motion occurs. 
This critical temperature is shifted to lower values for higher spin currents. In the limit 
of high effective spin current T* vanishes and domain wall motion can be observed over 
the whole temperature range. Here, the spin torque effect is no longer balanced by the 
anisotropy and only the terms responsible for the precession of the magnetization around 
the affects the domain wall motion. 

In both figures, analytical curves and numerical results agree. This demonstrates clearly 
that the assumption made for the derivation of Eq. [16] are reasonable for the parameters 
used. 

Note, that the pinning barrier which stops the domain wall motion can be overcome 
by thermal fluctuations. The role of these fluctuations was investigated by Duine et al.— 
within the framework of an extended, stochastic LLG equation. These fluctuations are 
relevant in the limit of very thin wires where by thermal activation the pinning potential 
can be overcome on sufficiently long time scales, leading to a finite domain wall motion even 
below the critical current. In our work, however, fluctuations are not considered so that the 
results are relevant for thicker wires where thermal fluctuations of the domain wall profile 
can be neglected. 
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FIG. 5: (Color online) The top figure shows the average domain wall velocity {v)g calculated from 
the LLB equation with Gilbert damping vs. the effective spin current u x for different temperatures. 
The bottom figure shows {v)q vs. the reduced temperature T/Tq. The lines represent the analytic 
solution (see Eq. [To]) and the points are from numerical simulations. The calculations are for the 
non-adiabatic case, (3q = A = 0.02. 

B. Non-adiabatic spin torque effects 



In the following, non-adiabatic spin torque is taken into account and its effect is discussed 
in more detail. For comparison with previous investigations,- 1 ^ the non-adiabatic prefactor 
Pg is assumed to be temperature independent and is investigated in relation to the tempera- 
ture independent microscopic damping constant A. Note, however, that in Ref. [3l| non-local 
contributions, which are strongly correlated to the domain wall width, are predicted, which 
are neglected for our wide walls here. 

Our results for the case (3q = A are shown in Fig. [51 In the zero temperature limit it is 
m = 1 and a± = A so that the last term of Eq. [10] vanishes and only the first spin torque 
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term remains finite which is responsible for the displacement of the domain wall along the 
x-axis. The magnetization is, hence, not tilted out of the easy plane and no torque occurs 
due to the precessional or relaxational part of the LLB equation. This behavior was already 
observed in previous numerical investigations of the LLG equation and is discussed in more 
detail in [(jl^. 

In the case of elevated temperatures, the situation is different due to the fact that the 
last term of Eq. [10] does not vanish because of the temperature dependence of m and a±. 
This term is responsible for tilting the magnetization out of the easy plane and it leads to 
the existence of the Walker threshold w Walkcr (Eq. fI5i) . Fig. [5] shows that two regimes can 
be distinguished: for u x < walker the velocity (v)q shows a linear behavior as in the zero 
temperature limit while in the regime u x > ^ Walkcr the second term in Eq. [16] takes over 
and the velocity increases even faster. Here, the last term of Eq. [10] leads to a continuous 
rotation of the magnetization around the x-axes. Following Eq. [15] the transition between 
these regimes is shifted to lower effective spin currents with increasing temperature vanishing 
at OK. 

The averaged domain wall velocity as a function of temperature for different effective 
spin current values is shown in the bottom part of Fig. [5] First, it decreases with increasing 
temperature until a minimum value is reached, after which the velocity increases. The 
minimum can be identified once again as the Walker threshold. This behavior is found for 
all effective velocities although the value of the minimum of the velocity is shifted to higher 
temperatures for lower effective velocities. 

Fig. [6] shows results for the case (3q < A. Here, the term responsible for tilting the magne- 
tization out of the easy plane plays a crucial role for all temperatures even at K. As before, 
the Walker threshold is shifted to lower effective velocities with increasing temperature. In 
comparison to the case (5q = A, the Walker threshold occurs at lower effective spin currents 
for the same temperature value so that the precession of the domain wall sets in earlier. 

The averaged domain wall velocity as a function of the temperature is shown in the 
bottom part of Fig. [6] for different values of the effective spin current. As before a minimum 
exists which is shifted to lower temperatures for higher effective velocities, consistent with 
the shift of the Walker threshold discussed above. In comparison to the (5q = A case, this 
shift of the minima is more pronounced. 

Finally, Fig. [7] shows results for the case (3q > A. Here, the behavior differs from the two 
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FIG. 6: (Color online) The top figure shows the average domain wall velocity {v)q calculated from 
the LLB equation with Gilbert damping vs. the effective spin current u x for different temperatures. 
The bottom figure shows {v)q vs. the reduced temperature T/Tq. The lines represent the analytic 
solution (see Eq. [T6|) and the points are from numerical simulations. The calculations are for the 
non-adiabatic case, /3q < A. 

cases before. First of all, above the Walker threshold, the wall velocity increases slower than 
linear not faster as before. This is due to the sign change in Eq. [161 Also, starting from low 
temperature the Walker threshold is first shifted to higher currents. In the effective spin 
current range shown in the figure the Walker threshold even disappears due to the fact that 
it is shifted out of the range presented. Surprisingly, at temperatures close to Tq another 
Walker threshold appears at lower effective velocities above which the averaged wall velocity 
increases faster than linear. 

This second transition can also be identified from the averaged velocity as a function of 
temperature as shown in the bottom part of Fig. [71 For an effective spin current of 600 m/s 
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FIG. 7: (Color online) The top figure shows the average domain wall velocity (v)g calculated from 
the LLB equation with Gilbert damping vs. the effective spin current u x for different temperatures. 
The bottom figure shows (v)q vs. the reduced temperature T/Tq. The lines represent the analytic 
solution (see Eq. ll6p and the points are from numerical simulations. For u x = 600 m/s two crossings 
of the Walker threshold are visible {T/Tq = 0.32 and T/Tq = 0.83). The calculations are for the 
non-adiabatic case, (3q > A. 

there are two temperatures where the behavior of the velocity changes. This corresponds 
to the fact that the Walker threshold as shown in Fig. [3] can be crossed twice for (5q > A 
and certain values of the effective current, leading to this intricate behavior, which could be 
easily identified if observed experimentally. 

C. Comparison of Gilbert and Landau-Lifshitz damping 

In this subsection the difference between the assumption of Gilbert damping on the one 
hand and Landau-Lifshitz damping on the other hand is discussed for the case (3q > A. As 
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FIG. 8: (Color online) The top figure shows the average domain wall velocity {v)g calculated from 
the extended LLB equation with Landau-Lifshitz damping vs. the effective spin current u x for 
different temperatures. The bottom figure shows vs - the reduced temperature T/Tq. The 

lines represent analytic solutions and the points are from numerical simulations. The calculations 
are for the non-adiabatic case, non-adiabatic case = 0.005 (/?g = 0.025) and A = 0.02). 

mentioned before the prefactors for the non-adiabatic spin torque term can be transformed 
as — PG~ a ±/ m c ■ This transformation is temperature dependent so that a qualitatively 
different temperature dependence exists for Landau-Lifshitz damping in comparison to the 
Gilbert damping discussed before. Fig. [8] summarizes our results for the extended LLB 
equation assuming Landau-Lifshitz damping (Eq. [S]). 

Comparing Fig. [7] and Fig. [8] it is found that the curves for K are indeed the same. 
However, the behavior at elevated temperatures is completely different: a Walker threshold 
can be identify for all temperatures. Below the threshold the domain wall is moved along 
the wire with an averaged velocity proportional to the effective spin current. Above the 
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threshold additionally to this movement the magnetization precesses around the x-axis, so 
that the velocity of the domain wall is decreasing. The temperature dependent threshold is 
shifted to lower effective velocities for higher temperatures. Here, the K curve is an upper 
limit for the averaged domain wall velocity. 

The behavior of the domain wall velocity as a function of the temperature is less compli- 
cated than in the corresponding Gilbert case. For low currents the averaged velocity steadily 
decreases with increasing temperatures. For slightly larger currents the domain wall moves 
with a constant velocity until the Walker threshold is reached. Above this point the velocity 
is decreasing. Close to the Curie temperature another transition point is reached where the 
domain wall velocity decreases even faster. Here, the Walker threshold is crossed again (see 
Fig. [3]). A further increase of the effective spin current leads to a shift of the first transition 
point to lower temperatures and only two regimes can be identified. 

VI. CONCLUSIONS 

In summary, we combined the LL form of the LLB equation as well as the Gilbert form of 
the LLB equation with the adiabatic and non-adiabatic spin torque terms. We investigated 
analytically as well as numerically domain wall motion at various temperatures for the 
adiabatic and non-adiabatic cases. The Walker threshold as well as the domain wall velocities 
show a strong temperature dependence. Furthermore, we found a different behavior for 
the temperature dependent Walker threshold assuming the Gilbert form of damping or 
the LL form. Since the two behaviors can be qualitatively different, a measurement of 
the temperature dependence of the velocity and the Walker threshold could pose a unique 
opportunity to identify, whether the Gilbert or the LL formulation are the physically relevant 
one. This then in turn could have implications for determining the physically relevant 
description of damping, which is one of the key open questions in magnetization dynamics. 
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